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<D ■ ABSTRACT 



Context. We study the precession of accretion disks in the context of gamma-ray burst inner engines. 
Aims. Our aim is to quantitatively estimate the characteristics of gravitational waves produced by the precession of the 
transient accretion disk in gamma-ray bursts. 

Methods. We evaluate the possible periods of disk precession caused by the Lense-Thirring effect using an accretion disk 
model that allows for neutrino cooling. Assuming jet ejection perpendicular to the disk plane and a typical intrinsic 
time-dependence for the burst, we find gamma-ray light curves that have a temporal microstructure similar to that 
observed in some reported events. The parameters obtained for the precession are then used to evaluate the production 
(■H ' of gravitational waves. 

Results. We find that the precession of accretion disks of outer radius smaller than 10 8 cm and accretion rates above 1 
I . Mq s _1 could be detected by Advanced LIGO if they occur at distances of less than 100 Mpc. 

Conclusions. We conclude that the precession of a neutrino-cooled accretion disk in long gamma-ray bursts can be 
probed by gravitational wave astronomy. Precession of the disks in short gamma-ray events is undetectable with the 
C/} , current technology. 
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■ 1. Introduction some that are difficult to explain with the standard model 

i^, _ . _ . /^.t^ \ ( e -g- Romero et al. 1999). It has been suggested that the 

l In the central engines of gamma-ray bursts (GRBs), ac- precession of the jet can pl a role in the formation of the 

^ . cretion onto a black hole resulting from the collapse of a microstructure of both long and short gamma . ray bursts 

. ■ massive star (e.g. Woosly 1993) or a merger of two com- ( Blackman et al. 1996, Portegies et al. 1999, Fargion 

w a | pact objects (e.g. Mocnkovitcn et al. 1993) leads to the lor- -^ggg R e y noso e ^ a i 2006) 

, mation of a hot and dense, transient accretion disk. This ' , , i i i ^ ^ r 

O-,. n i . -n iii j. • i rr^i Reynoso et al. (2006) developed a model lor precessing 

disk can be significantly cooled by neutrino losses. I he ac- . , , J , . . v , y , . , ° 

,. c i, , , m 7\/f -l • lets based on spin-induced precession of the neutrino-cooled 

cretion of matter, at a rate ~ 0.1 — 10 M ft s , is sup- J . . ^ . f , . J ... , , 

j , -i i , j xi j • . • . ,i massive disk. Ine precession ot the disk is transmitted to 

posed to power the burst, and the radiative processes m the , ,. . ,. . ^ ^ u . . ±1 , 

n ,. i i i j. r u j the relativistic lets, resulting m the peculiar temporal mi- 

l j relativistic lets are expected to account for the observed J J / ^t-.t-» 5 t i , 

. v , , r^i v i -1 r crostructure ot some GRB s light curves. 

rN light curves. Ihese curves display a wide variety ol time e it rc 

b profiles with timescales from milliseconds to minutes. The In thl r s P a P er ' we stud ^ an additional effect of this pre- 

■ - - usual interpretation of this temporal structure is in terms of cession of ver ^ massiv ^ accretion disks: the production of 

shocks that convert bulk kinetic energy into internal energy gravitational waves. Gravitational wave radiation is ex- 

of the particles, which then cool by means of synchrotron P ected fr ? m gamma-ray bursts if the gravitational collapse 

and inverse Compton emission. The shocks can be internal is non-spherical, if there are strong inhomogeneities in the 

to the jet and produced by colliding shells with different accretion disk, or, in the case of short GRBs as the re- 

Lorentz factors (e.g. Kobayashi, Piran & Sari 1997; Daigne sult ^ the spiraling and merging of compact objects (e.g. 

& Mochkovitch 1998; Guetta, Spada & Waxman 2001) or Mineshinge et al. 2002, Kobayashi & Meszaros 2003). 
the result of interactions with the ambient medium (e.g. The S1 g nal we consider here has a different origin from 

Heinz & Begelman 1999; Fenimore, Madras & Nayakshin those previously discussed in the literature, and its specific 

1996). Among the observed light curves, however, there are features can shed light on the behavior of the innermost 

regions of the sources. We demonstrate that gravitational 

Send offprint requests to: G. E. Romero wave astronomy with instruments such as Advanced LIGO 

e-mail: romero@iar-conicet .gov. ar can be used to probe the Lense-Thirring effect in nearby 
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2. Accretion disks and spin-spin interaction in GRB 
engines 

Transient accretion disks are formed in GRB's engines such 
as collapsars and mergers of compact objects. The accretion 
rate in these disks is expected to vary significantly mostly 
in the outer part of the disk, while for the inner disk a 
constant accretion rate remains a valid approximation (e.g. 
Popham 1999, Di Matteo et al. 2005). 

The cons ervation of mass falling with a velocity v r ~ 
ry/GM^r -3 at a radius r from the black hole axis is given 

by 



M = -27rrw r E(r), 



(1) 



where E(r) = 2p(r)H(r) is the surface density, H(r) is the 
disk half- thickness, and p(r) is the mass density of the disk. 
The conservation of angular momentum and energy can be 
used to obtain numerically the functions S(r) and H(r), 
assuming that the heat generated by friction can be bal- 
anced by advection and neutrino emission (Reynoso et al. 
2006). If the orbit of a particle around a spinning black 
hole is not aligned with the black hole equator, then the 
orbit precesses around the spinning axis. This is called the 
Lense-Thirring effect (Lense & Thirring 1918), and origi- 
nates from the dragging experienced by the inertial frames 
close to the rotating black hole. In accretion disks, the ac- 
tion of viscous torques leads to the alignment of the very 
inner part of the disk with the black hole equator plane 
(Bardeen & Petterson 1975). 

If the transient accretion disk in GRBs is formed mis- 
aligned with respect to the equator of a rapidly spinning 
black hole, then the accretion disk will develop precession. 
It has been argued that the disk precesses approximately 
like a rigid body, i.e., it does not present warping, when the 
disk Mach number is M. < 5 (Nelson & Papaloizou 2000). 
This condition is fulfilled in the accretion disks of GRBs. 
The precession of the disk leads to the precession of the jets, 
yielding a source of temporal variability and microstructure 
in the signal (Reynoso et al. 2006). 

Neglecting any nutation movement, the precession pe- 
riod of the disk r p can be related to its surface density as 
(see Liu & Melia 2002, Caproni et al. 2004, Reynoso et al. 
2006): 



2tt 



d sin0d(j) = 27rsin#7 d 
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where the magnitudes of the disk angular momentum 
and the precessional torque applied to the disk are 

rRout 



= 2tt [ 
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Here, if a* is the spin parameter, 
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is the relativistic Keplerian angular velocity, R g = 
GMbh/c 2 is the gravitational radius, and 



2tt 



1 
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(6) 



is the nodal frequency obtained by perturbing a circular 
orbit in the Kerr metric (Kato 1990). The precessing part 
of the disk ends at an outer radius R ou t , extending from an 
inner radius R ms = £ ms ^g, where 



6ns = 3 + A 2 t [(3 - Ai)(3 + A 1 + 2A 2 )] X / 2 , 
with 

A x = 1 + (1 - alf/\l + a,) 1 / 3 + (1 - a,) 1 / 3 ] 
and 

A 2 = {Zal + A\y'\ 



(7) 
(8) 
(9) 



The minus sign in the expression for £ ms corresponds to pro- 
grade motion (a* > 0), whereas the plus sign corresponds 
to retrograde motion (a* < 0). 

In Fig. [TJ we show the precession period obtained as 
a function of the disk outer radius for different accretion 
rates. 
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Fig. 1. Precession period as a function of the outer radius 
of the accretion disk for M = O.1M s _1 (dotted line), 
M = 1M© s -1 (solid line), and M = 10M o s -1 (dashed 
line). 



3. Gravitational waves from precessing disks in 
GRBs 

We attempted to characterize the emission of gravitational 
waves (GWs) from precessing accretion disks in GRB en- 
gines. Precession leads to observable features in the gamma- 
ray light curves. Its effects are illustrated in the following 
subsection, and the necessary formulae for GW emission 
are presented in Sect. 3.2. 

3.1. Specific models of precessing accretion disk in GRBs 

We assumed that the precession of accretion disks in GRBs 
is transmitted to the jets. The electromagnetic emission is 
generated in the jet by synchrotron and inverse Compton 
processes. The dependence of the gamma-ray luminosity on 
the angle with respect to the jet axis ip is described by the 
expression (Portegies-Swart et al. 1999) 
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L(V) 



27 

4 



-0.6*(x) 



d^(x) 



(10) 



-0.3^0) _ e -6.3^(cc) 

J dx 

where x = 10 sin^ and ^(x) = | ln(l + 4x 2 ). The intrinsic 
temporal dependence of the signal is characterized by a 
FRED (Fast Rise and Exponential Decay) behavior, 



/to 



tan 



t 



7"plat 



(ii) 



where Ni is a normalization constant such that the max- 
imum of the signal corresponds to unity, and r r i se , T p i at , 
and Tdec are the timescales for the initial rise, plateau, and 
decay, respectively. 

We considered two specific GRBs as examples: the 
short burst GRB 990720 and the long burst GRB 990712. 
We found for these bursts, a proper set of the relevant 
timescales and precession period. Then, using the model 
of Reynoso et al. (2006) we reproduced the observed light 



Table 1. Parameters used for the example GRB events. 



Parameter GRB 990712 GRB 990720 



Tprec [s] 
Trise [s] 
Tplat [s] 
Tdec [s] 

M[M©s _1 ] 
i?out[10 7 cm] 



1.0 
2.0 
10.0 
4.8 
1 - 10 
3.3-3.7 



0.2 
0.04 
0.5 
0.14 
0.1 
2.1 



The original time profiles corresponding to the men- 
tioned events and the light curves obtained using F(t) = 
I(t)L(ij)(t)) where the effect of precession is given by (j)(t) = 
27r(t/r pre c), are shown in Figs. [2] and [3l In both cases, the 
observer is located at o bs = 2° with respect to the z-axis 
perpendicular to the black hole equator, and the angle %jj{t) 
between the jet and the observer is time-dependent because 
of the precession. 



BATS E Trigger 7653 




Since Trigger (9907Z0 : 75940.545) 




Fig. 2. GRB 990720 light curve (left panel) and light curve 
obtained with a precessing disk (right panel). 
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Fig. 3. GRB 990712 light curve (left panel) and light curve 
obtained with a precessing disk (right panel). 

3.2. General formulae for gravitational wave emission 

An axissymmetric body (i.e., with inertial moments I\ = 
I2) in precession emits gravitational waves with an ampli- 
tude given by (Zimmermann & Szedenits 1979, Maggiore 
2008) 



Krec(t) = h+(t) + h x (t), 

where 

h+(t) = cos Qt H 
h x (t) = F Xj i sin lit - 

with 



F+, 2 cos2fit, 
- F x 2 sin 2^t, 



and 

K = 



= h f sin 2a sin i cos i 



(12) 

(13) 
(14) 

(15) 
(16) 
(17) 
(18) 

d ' (19) 
Here, a is the angle between the angular momentum of the 
disk and that of the black hole, i is the angle between the 
z-axis of the detector and the signal direction of arrival, 
and d is the distance to the radiating body. The principal 
inertia moments are 



2/ig sin 2 a(l + cos 2 l) 
h f sin 2a sin i 
4/ig sin 2 a cos l. 



G(i 3 -h)n 2 



(x 2 



h 



-I 



■ y 2 )p{r) 
y 2 )p(r)- 



(20) 
(21) 



The frequency of the gravitational waves are f\ = fi/(27r) 
and /2 = 2^/(27r), which are related to the angular mo- 
mentum of the body by 



(22) 



Since we are dealing with bursting events, the GW sig- 
nal is expected to be significant for a time T p i at . Therefore, 
we modulated the signal of Eq. ([T2]) using a Gaussian 



h(t) = h prec (t) e 



7?— 

plat 



(23) 
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This is usually adopted to describe the GW signal from 
bursting sources (e.g. Acernese 2008, Maggiore 2004). The 
angular frequencies that contribute to the waveform of Eq. 
([23]) can be obtained from its Fourier transform h(uj) = 
h+(uj) + h x (cj), where 



and 

h x (uj) 



Tplat 



F. 



+,2 e 



plat 



Tplat 



F x i ( e^ 1 ^ 2 



-|^(w-2n) 



T plat 



(24) 



• e 2 



(uj-ny 



Fx,2 (e 



. (25) 



It can be seen from these expressions that the ranges of fre- 
quencies that contribute to the signal are centered around 
Q and 20, and the width of the ranges is r~^ t . Given the 
typical durations of GRBs, the frequency spread is narrow 
for all bursts except for those with durations much shorter 
than 1 s. In Fig. H]we show an example waveform obtained 
with Eq. ([23]) for an accretion rate M = lM s _1 , a pre- 
cession period r prec = 0.3 s, r p i at = 10 s, i = 45°, a = 20°, 
and d = 100 Mpc. 
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Fig. 4. A typical GW waveform that can be produced by a 
precessing GRB accretion disk. The accretion rate is M — 
IMqS -1 , the precession period is T prec = 0.3 s, and the 
duration of the event is given by the timescale r p i at = 10 s. 



4. Detectability 

We now consider the gravitational waves produced by the 
disk precession in events with the characteristics of those 
discussed in the previous section. In Table 1, we list the 



main parameters for the short and long GRBs taken as 
examples. We estimated the amplitude of the gravitational 
waves using Eq. ([T9]h considering different accretion rates. 
In order to assess the detectability, we calculated the root- 
sum-square amplitude (e.g. Acernese 2008, Maggiore 2004) 



dt{h\ + h 2 x ), 



(26) 



where / = Q/(2ir). For illustration, we choose the value 
l = 45° for the angle between the line of sight and the z-axis 
of the detector, and we take a = 10° — 20° to be the angle 
between the angular momentum of the disk and that of the 
black hole. In Fig. [5] we plot the root-sum-square amplitude 
h T88 , as a function of one of the allowed frequencies, / = 
17/ (2tt), and also as a function of the outer radius of the 
inner precessing disk. In the figure, we include the expected 
sensitivity for Advanced LIGO (Shoemaker 2010). 

The parameters used to find E(r) and H(r) are 
Mbh = 3M , a* = 0.1, a viscosity parameter 
a = 0.1, and the different mass- loss rates M = 
{O.1M s- 1 ,lM s- 1 ,10M o s -1 } . The distance to the 
GRB is taken to be d = 100 Mpc. 

As can be seen from Fig.[5j there are higher probabilities 
of detection for accretion rates higher than 1M s _1 and 
outer radii between 10 7 and 10 8 cm. When the accretion 
rate is very high, the disk may become advection dominated 
rather than cooled by neutrino emission (Liu et al. 2008). 
The dynamics in the gravitational field, however, is not 
affected. High accretion rates can be sustained only in long 
GRBs, so we conclude that there is only a good prospect 
for detection of gravitational waves from precessing disks 
of nearby (d < 100 Mpc) and long events. These events are 
likely to be related to the death of very massive stars, so 
the host galaxies should have active star- forming regions. 

Two low-luminosity long GRBs (980425 and 060218) 
were already observed at distances of ~ 40 Mpc and ~ 
130 Mpc, respectively (Corsi & Meszaros 2009). The local 
rate of long GRBs is estimated to be ~ 200 Gpc -3 yr _1 
(e.g. Liang et al. 2007, Virgili et al. 2009). INTEGRAL has 
detected a large proportion of faint GRBs that have been 
inferred to be local (Foley et al. 2008). All this suggests 
that the detection of precessing disks of GRBs through their 
gravitational emission is possible in the near future. 

5. Summary and conclusions 

We have studied long and short gamma-ray bursts showing 
that the presence of microvariability in their light curves 
can be caused by the precession of the transient accretion 
disk. We then considered precessing disks of GRBs as a 
source of GWs and estimated the resulting waveform pro- 
duced by the phenomenon. 

Our results indicate that if the outer radius of the pre- 
cessing disk is between 10 T and 10 8 cm and the accretion 
rate higher than 1 M s _1 , the gravitational waves can be 
detected from distances of 100 Mpc or less. We conclude 
that only in relatively nearby and long GRBs can the pre- 
cession of the accretion disk be studied by gravitational 
wave astronomy. The detection of one event of this class 
can be used to test the Lense-Thirring effect in the strong 
field limit. 
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Fig. 5. Gravitational wave rss amplitude for different accretion rates and Advanced LIGO sensitivity (red dashed line) as 
a function of the gravitational wave frequency for a = 10° (left panel) and for a = 20° (right panel). The corresponding 
outer radius of the accretion disk is indicated in the upper horizontal axis. 
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